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Abstract 

We give a general class of static, spherically symmetric, non-asymptotically flat and asymp- 
totically non-(anti) de Sitter black hole solutions in Einstein-Maxwell-Dilaton (EMD) theory of 
gravity in 4-dimensions. In this general study we couple a magnetic Maxwell field with a general 
dilaton potential, while double Liouville-type potentials are coupled with the gravity. We show 
that the dilatonic parameters play the key role in switching between the Bertotti-Robinson and 
Reissner-Nordstrom spacetimes. We study the stability of such black holes under a linear radial 
perturbation, and in this sense we find exceptional cases that the EMD black holes are unstable. 
In continuation we give a detailed study of the spin-weighted harmonics in dilatonic Hawking ra- 
diation spectrum and compare our results with the previously known ones. Finally, we investigate 
the status of resulting naked singularities of our general solution when probed with quantum test 
particles. 
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I. INTRODUCTION 



We revisit the 4— dimensional Einstein-Maxwell-Dilaton (EMD) theory and show that 
there are still plenty of rooms available to contribute the subject. Double Liouville potential 
and general dilaton coupling is considered to obtain more general solutions with extra pa- 
rameters and diagonal metric in the theory. From the outset we remind that, depending on 
the relative parameters, the double Liouville potential has the advantage of admitting local 
extrema and critical points. The Higgs potential also shares such features, whereas single 
Liouville potential lacks these properties. Double Liouville-type potentials arise also when 
higher-dimensional theories are compactified to 4— dimensional spacetimes and expectedly 
bring in further richness. All known solutions to date can be obtained [l-3] as particular 
limits of our general solution, and it contains new solutions as well. In the most general form 
our solution covers Reissner-Nordstrom (RN) type black holes and Bertotti-Robinson (BR) 
spacetimes interpolated within the same metric. Interpolation of two different solutions in 
general relativity is not a new idea J4]. Particular limits of the dilatonic parameter yield 
the RN and BR spacetimes. In between the two, the linear dilaton black hole (LDBH) lies 
for the specific choice of the parameters. It is well-known that the near horizon geometry 
of the extremal RN black hole yields the BR electromagnetic universe. The latter js| is 
important for various reasons: It is a singularity free non-black hole solution which admits 
maximal symmetry and finds application in conformal field theory correspondence (i.e. AdS/ 
CFT). Particles in the BR universe move with uniform acceleration in a conformally flat 
background. These features are mostly valid not only in N = 4 but in higher dimensions 
(N > 4) as well. The topological structure of the BR spacetime is still AdS2 x S N ~ 2 in 
N— dimensions with the radius of S N ~ 2 depending on the dimension of the space. Recently 
we have extended the Maxwell part of the BR spacetime to cover the Yang-Mills (YM) field 
and obtained common features that share with the Maxwell field {(J. The dilatonic black 
hole solution involved in the general solution obtained in this paper is non-asymptotically 
flat, therefore we expressed it in terms of the quasi local mass (Mql) Q. The metric is reg- 
ular at horizons with only available singularity at r = 0. Another feature is the asymptotic 
(r — > 00) absence of (anti) de-Sitter property which was discovered also within the context 
of different models 3']. Our general solution has been tested for stability against the radial, 
linear perturbations. We found that presence of dilaton can trigger instability in the RN 
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black hole which is stable otherwise. Our analysis proves that the BR sector remains man- 
ifestly stable against such perturbations. Thermodynamic stability has also been discussed 
briefly by considering the specific heat of the metric. Divergence in the specific heat for 
specific values of the parameters signals phase transition in our thermodynamic system, i.e., 
topology change in the spacetime. 

Next, we concentrate ourselves on the LDBH case and analyze the Hawking temperature 
both from semi-classical and standard surface gravity methods We point out the 

contrasts between the two methods when there are single and double horizons. The high 
frequency limit of the semi-classical radiation spectrum method (SCRSM) does not agree 
with the Hawking's result. It is observed, as an interesting contribution in this work that 
the coupling between scalar field charge and the magnetic charge of the spacetime gives rise 
to spin- weighted spheroidal harmonics which plays a dominant role in the difference. In the 
absence of such coupling, when the scalar field is assumed chargeless for instance, similar 
analysis was carried out previously and we had recovered the same results easily. It turns 
out that the very existence of a spin-weighted spheroidal harmonics in the theory transforms 
a divergent temperature spacetime to a finite one. We argue that such a behavior may play 
a leading role in the detection of such LDBHs. 

In the final section of the paper we appeal once more to the test scalar field equation, but 
this time with the purpose to investigate the quantum nature of the naked singularities. We 
identify first the particular solution that yields horizonless naked singularity at r = 0. By 
invoking the Horowitz-Marolf 10( criterion on quantum nature of classical singularities we 
explore under which set of parameters classically singular but quantum mechanically regular 
metrics can occur in our general solution. 

The organization of our paper is as follows. In Sec. II we introduce our action, field 
equations and obtain the general solution. Sec. Ill singles out the linear dilaton case 
and investigates the stability of our general solution. Application of the SCRSM and its 
connection with the Hawking temperature is employed in Sec. IV. Sec. V discusses the 
status of naked singularities from quantum picture. We summarize our results in conclusion 
which appears in Sec. VI. 
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II. FIELD EQUATIONS AND THE METRIC ANSATZ FOR EMD GRAVITY 



The 4— dimensional action in the EMD theory is given by (87rG = 1) 

S = J d 4 x^ Qi? - l -d,4>d^ - V (0) - l -W (0) (F Xa F x ^ , (1) 

where 

V ((f)) = V x e M + V 2 e M , W (0) = A^ 271 * + A 2 e~ 2720 . (2) 

refers to the dilaton scalar potential and 73 denotes the dilaton parameter, Aj is a constant 
and V (0) is a double Liouville-type potential. We note that we exclude the simultaneous 
values j3i = 02 and 71 = 72 in general, since these particular values lead to the already 
known cases. 

Let us remark that although double Liouville potential in V (0) , which renders local 
minima, necessary for construction of vacuum states possible, the similar choice for W (0) 
seems less appealing. It will be justified from the exact solutions below, however, that there 
are asymptotics which remains inaccessible by the choice of a single Liouville term in W (0) . 
Stated otherwise, at both asymptotes of r = and r = 00 (or f = and f = 00 for LDBH) 
dilatonic coupling to the magnetic field becomes much stronger. Choosing a single Liouville 
potential simply looses the strength at one end of the range. Besides, it is all a matter of 
choice to set Ai (A2) = 0, which makes the dilatonic coupling asymptotically free. In the 
LDBH case as it will be proved, if we set Ai = 0, we shall remove the possibility of an inner 
(Cauchy) horizon which justifies the advantages and motivation for choosing the double 
Liouville-type potential in W (0) . In (1) R is the usual Ricci scalar and F = ^F^dx^ A dx v 
is the Maxwell 2— form (with A indicating the wedge product) given by 

F = dA, (3) 

for A = A^dx^, the potential 1— form. Our pure magnetic potential with charge Q, which 
is given by 

A = -QcosO d(p, (4) 

leads to 

F = Qsm8 d6Ad(p. (5) 

Let us note that with the present choice of W (0) the electric-magnetic symmetry that exists 
in the standard dilatonic coupling, i.e., Ai (A 2 ) = 0, is no more valid. Our choice in this 
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paper relies entirely on the magnetic choice. Variations of the action with respect to the 
gravitational field and the scalar field lead, respectively to the EMD field equations 



R,u = d M <R0 + V (0) g, u + W (0) (2F^F V X - l -F Xa F x °g,)j , 



(6) 



V 2 - V (0) - -W (0) (F Xa F x °) = 0, (7) 

_ d\ 
/= d0y!' 

where is the Ricci tensor. Variation with respect to the gauge potential A yields the 
Maxwell equation 



d(W(0)*F) = O, (8) 
in which the hodge star * means duality. 

A. Ansatz and the Solutions: 



Our ansatz line element for EMD gravity is chosen to be 

ds 2 = -f (r) dt 2 + -r^-rdr 2 + R(r) 2 (d6 2 + sin 2 6dip 2 ) , (9) 
with / (r) and R (r) only function of r while the Maxwell invariant takes the form 

F»F*" = ^. (10) 
The Maxwell equation (8) is satisfied automatically and the field equations become 

vV : = (fl 2 yy)' = V (</>) + \w W (f^), (n) 

fl ^j^^ w _E^W n (12) 

« ;:= _^_(^ =/ ^ w _m (FA ^), (13) 



in which a prime stands for derivative with respect to the argument of the function. We 
start with an ansatz for R(r) as 

R(r)=Ae^ (15) 
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in which A and r\ are constants to be found. Substitution in (12) and (13), implies 

2r] 



0(r) 



2r] 2 + l 



In r. 



(16) 



Finally by putting these results into Eq.s (11) and (14) one finds that by setting 

1 



V 



a 



V2' 



(17) 



and 



a 



1 



71 = y/2' 72 ay/2' 

1 

a 



(18) 



ft = V2a, fa 



a general solution for /(r) reads 

/(r) = (l + a 2 ) 2 



g2 AirT ^ 
(1 + a 2 ) A 4 



+ 



(19) 



Q 2 A 2 
A 4 



with the constraint condition 



-y 2 



7" !+° 



'1 + a 2 ) a 2 3 — « 2 



— Mr i+" 2 



V 2 (1 - « 2 ) A 4 - a 2 A 2 + A 2 g 2 (1 + a 2 ) = 0. 



(20) 



Herein M is a mass-related integration constant and a and A are constants that will serve 
to parametrize the solution. We note that in case that we are interested in the Newtonian 
limit, when a = = X 2 = V 2 , and r — > oo, we must choose M — > 2M, so that M represents 
the Newtonian mass. Therefore the dilatonic function 0, Liouville potential V and in 
terms of a become 

ay/2 



<f>(r) 



1 + a 2 



lnr, R (r) = Ar 1 + a ' 2 , 



\/ = \/ir + V 2 r , = Airi+^ + A 2 n+« 2 . 

We remark that the solution (16-20) is the general diagonal solution that covers all particular 
solutions of this kind known so far. For arbitrary value of a, other then 0, 1 and oo, it yields 
a new solution in accordance with our ansatz. It is observed also that our metric and 
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potentials are invariant under a — > —a, whereas — > — 0. The asymptotic behavior of the 
metric function, / (r) and other limiting cases can be summarized as follows 



2\ 2 



(l + a 2 ) 



lim / (r) — > < 



y ir l+o^ 
3-a 2 



2 2 A 2 



-Fi-y 2 r 



< a 2 < 1 
2 -l 

2 



(21) 



[{&) (^-V 2 )r^Ka 
lim / (r) ->■ (1 + a 2 ) ( 



Q 2 Xi 



(22) 



The case a 2 = 1 will be studied separately, while the case a 2 = 0, with the choice of A 2 = 0, 
V 2 = 0, and A — 1 — X 1 leads to 



/(r) = l-— r - — + -2-, 



i?(r) = r,0 = 0, 



(23) 



which corresponds to the action 



(24) 



This is recognized as the 4— dimensional action in the EM theory with the solution rep- 
resenting a RN black hole with a cosmo logical constant. Another limiting case of interest 
consists of the case with a 2 — > oo, A 2 = 1, with the action 



(25) 



leading to the solution 



(26) 



f(r)= r 2 -Mr, 

A 2 (V 2 A 2 + 1) = Q 2 , 
R(r)=A, 0(r) = O, 



in which M is the mass related integration constant. Here also we have a 4— dimensional 
action in the EM theory with cosmological constant but the metric function represents a 
BR space time. 

By looking at the asymptotic behaviors of the general solution one finds that < a 2 < 1 
and 1 < a 2 correspond to the cases of RN and BR solutions, respectively. Here a 2 = 1 acts 



much like a phase transition which changes the structure of space time from RN into BR. 
The thermodynamic instability from the expression of specific heat capacity Cq, (Eq. (58) 
given below) justifies this fact. It is quite interesting to see what will be the answer if one 
chooses a 2 = 1. In the next section we concentrate on this critical value for a 2 . 



III. THE LINEAR DILATON 

From the asymptotic behavior of the metric function one may see that a 2 = 1 is a critical 
value and the behavior of spacetime changes. In this chapter we only concentrate on this 
specific value for a 2 , and will be referred to as linear dilaton. The general solution after this 
setting reads 

"71=72 = ^, (27) 
ft = (3 2 = V2, 

4>(r) = — /=lnr, R(r) = Ay/r, 
v2 

V Ai 
V=-, W = ^ + X 2 r, 

A 2 = 2\ 2 Q 2 , (A 2 > 0) 



f(r) = 



A +(^-2V) r-M 



\ 2 A 2 r \A 2 



(28) 



where V = V 1 + V 2 . 

In order to explore the physical properties of the linear dilaton case we perform the 
transformation R (r) = A\/r — > r. This transforms the metric into, 

ds 2 = -f{¥)dt 2 + -j^d^ 2 + ?W, (29) 

in which 

/w4((i-^)5- M ^ + |)' < 30 » 

where the mass Mql denotes the quasilocal mass whose general definition is given below in 
Eq. (46). Other related parameters transform into the following forms, 
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4,(7) = V2]n(j\, (31) 
V(7) = * (Vl + V2) , W + 

The location of horizons can be found if we set the metric function g tt = 0. The solution is 

r h = ^=^M QL ±^M 2 QL -Aac, (32) 

where 

The linear dilaton solution admits single or double-horizons if the parameters are chosen 
appropriately. Another possible case is the extremal limit that occurs if Mq L = Aac. The 
horizon in this particular case is given by = y ^n^- The double horizon case occurs if the 



parameters simultaneously satisfy Mql > Mq L — Aac and Mq L > 4ac. This choice leads 
to the horizons 



\ 



Mql + ^ 



Ml 



QL 



Aac 



2a 



(34) 



M, 



QL 



M 2 QL - 



2a 



The relations between the parameters and double-Liouville-type potentials in the forma- 
tion of black holes becomes evident if one looks for the critical case. This is the case when 



Mql = yj Mq L — Aac, which follows that 2V = Hence if 2V < ^j, no horizon forms 
and the central singularity r = becomes a naked singularity. It can easily be seen that for 
Ai = 0, or for the single Liouville-type potential in W ((f)) , we have automatically single, 
outer event horizon alone. Another interesting property is in the behavior of the curvature 
scalar R. The curvature scalar for the metric function (29) is, 



R 



4f 4 (aA 4 - 1) - A A (or 4 - Mqlt 2 + c) 
2?e 



(35) 
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Note that the curvature scalar is finite at the location of horizons. Furthermore, when 
r — > oo , the Kretschmann and curvature scalars, the Liouville-type potentials and the 
coupling term of dilaton with Maxwell field all vanish. The mass and charge are finite and 
the dominant field is gravity with finite curvature. Consequently, the solution given in Eq. 
(29) is well-behaved. However, the Q = limit does not exist. 



A. Linear stability analysis of the general solution 



By employing a similar method used by Yazadjiev [11] we investigate the stability of the 
possible EMD solution, in terms of a linear, radial perturbation. To do so we assume that 
our dilatonic scalar field <fi (r) changes into O (r) + if) (t, r) , in which if) (t, r) is very weak 
compared to the original dilaton field <p (r) and we call it the perturbed term. As a result 
we choose our perturbed metric as 

ds 2 = -f (r) e r ^dt 2 + e x( *' r) -^y + R (rf (M 2 2 . (36) 

One should notice that, since our gauge potentials are magnetic, the Maxwell equations 
(Eq.(8)) are satisfied. The linearized version of the field equations (11-14) plus one extra 
term for R tr are given by 

Rtr : Xt = d r <P (r) d t if> (t, r) (37) 

- X V 2 J + ~ (T - x) r tif - 91 V (0 O ) ip = -^d 2 o W (0 O ) V (38) 
* R(r) 

R ee : (1 - Rose) X ~ \RR 7 (r - x) r = (r% V (<f> ) + j^d^W (<f> )) if; (39) 

in which a lower index represents the quantity in the unperturbed metric. First equation 
in this set implies 

X{t,r) = -t/){t,r) (40) 
V 

which after making substitutions in the two latter equations and eliminating the (r — x) r 
one finds 

V 2 ip(t,r) -U (r)iP(t,r) =0 (41) 

where 
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To get these results we have implicitly used the constraint (20) on A. Again by imposing 
the same constraint , one can show that U (r) is positive. It is not difficult to apply the 
separation method on (41) to get 

*P (t, r) = e^C (r) , V^C (r) - U eff (r) C (r) = 0, U eff (r) = (j + U (r)Y (43) 

where e is a constant. If one shows that the effective potential U e ff (r) is positive for any 
real value for e it means that there exists a solution for £ (r) which is not bounded. In other 
words by the linear perturbation our black hole solution is stable for any value of e. 

But in our case one must be careful. For instance let's go back to the general solution 
(19) and set V x = 0, 

this solution may have double horizons, single horizon (extremal) or no horizon. These 
depend on the values of the parameters. One may notice that this solution is a non- 
asymptotically flat metric and therefore the ADM mass is not defined in general. Following 

n 

the quasilocal mass formalism introduced by Brown and York |7j it is known that, a spher- 
ically symmetric N— dimensional metric solution as 

dR 2 

ds 2 = -F {Rf dt 2 + ——5 + R 2 dQ 2 N _ 2 , (45) 

G [H) 



aa 



admits a quasilocal mass Mq L defined by 

N — 2 

M QL = — g— <- 3 ^ (Rb) (G ref (Rb) ~ G (R B )) . (46) 

Here G re f (R) is an arbitrary non-negative reference function, which yields the zero of the 
energy for the background spacetime, and Rb is the radius of the spacelike hypersurface 
boundary. Applying this formalism to the solution (44), one obtains the horizon M in terms 
of M QL as 

M = jrrh)^ M ^ (47) 

after which the metric function becomes 

Indeed, since we wish to cover all known solutions in the literature of this kind, we consider 
Ai, M QL > 0, and 



Q 2 \ 



A 4 
11 



2 - V 2 ) > 0. (49) 



This condition together with Eq. (20) give a transparent view of U e ff (r) . In other words, 
after simplification, one can rewrite U (r) as 

which reveals for — < V 2 < U (r) and then U e ff (r) are positive, which means 

that the corresponding metric is stable. But for V 2 < — ^jt 2 , if ct 2 < o? critical where 



I Vol - 



2 A 2 



2 I ' z I A 4 / cr i \ 

^critical , T ~ | o 2 Ao ' V / 

|F 2 | + ^ 

then U (r) gets negative value and therefore our solution faces an instability condition. Here 
it is interesting to note that a 2 critical < 1 belongs to the RN type black hole solutions, i.e. 
BR type solution is automatically stable for any value of a 2 . 

The general solution reveals another interesting case after we set V\ = 0, and A 2 = i.e. 

, w -fi^{-^-^ + <£}. 

Upon choosing V 2 < (V 2 > 0) this admits the effective potential 

U{r) M^l) ,53) 



which clearly from (20), for a 2 < 1 (a 2 > 1) manifests an unstable black hole solution. As a 
result we observe that a stable RN black hole becomes unstable under certain conditions in 
the presence of a dilaton and a Liouville potential. 



B. Thermodynamic stability 

Concerning the solution (19), we set the parameters Ai = A 2 = 1 and VI = to get 
,. (l+a 2 ) \(Q 1 \ r 2 2M QL Q 2 \ 

which in terms of the radius of horizon r h one finds the quasilocal mass as 

_ r\ (Q 2 - V 2 A*) + QW 
Mql ~ 2AWF h • (55) 



The Hawking temperature 

/' (r h ) _ (1 + a 2 ) [r 2 (A 2 - 2Q 2 ) - Q 2 (1 - a 2 )] 



■ H 



47T 2±4 

A(l-a 2 )AHrf^ 
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(56) 



and the Bekenstein-Hawking entropy 



a 

4 



S = - = nr 



(57) 



where a is the area of the black hole, together lead to the heat capacity Cq for constant Q 
as 



Cq — T H 



dS 



dT, 



(a 2 + 1) 









+ l-a 2 


(a 2 - 1) 








+ 3 + a 2 



2nri 



(58) 



Our black hole solution becomes thermodynamically stable /unstable depending on Cq > 
/Cq < which is not difficult to test from this expression. For C^j < 2 and a 2 < 1, as an 



example, our black hole becomes thermodynamically unstable. Also for 



2 one gets 



Cq 



a 2 + l 



2nr 



3 + a 2 

which shows an instability independent of the values of a. Tab. 1 illustrates the stable and 

^ ^ ■» ,™ * (2 -($)')• 





a 2 — 1 < x 


-1 < x < a 2 - 1 


- (3 + a 2 ) < x < -1 


x < - (3 + a 2 ) 


a 2 < 1 


Unstable 


Stable 


Stable 


Unstable 


a 2 > 1 


Stable 


Unstable 


Unstable 


Stable 



(Table: 1) 



Eq. (58) reveals also that a 2 = 1 (i.e., the linear dilaton) is a phase transition point, 
however, there may be other possible transition points following a solution for a in the 
quadratic equation 



rl\2 



+ 3 + a 2 = 0. 



(59) 



IV. APPLICATION OF THE SCRSM AND HAWKING TEMPERATURE 



In this section, we shall attempt to make a more precise temperature calculation for 
the non-extreme LDBHs, a 2 = 1 given in Eq. (28), by using a method of semi-classical 
radiation spectrum, which has been recently designated as SCRSM 9] . The main difference 
between our present work with others 0, 9] (and references therein) is that the considered 
non-extreme LDBHs possess two horizons, due to having magnetic charge, instead of one. 
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Here, we first consider a massless scalar field \I/ with charge q obeying the covariant 
Klein-Gordon equation in the LDBH geometry. Namely, we look for the exact solution of 
the following equation, 



□* = 0, 

where the d' Alembertian operator □ is given by 



(60) 



□ 



1 



-g 



(61) 



in which symbolizes the covariant gauge differential operator as being 



Dp = — iqA„. (62) 

The scalar wave function \1/ of Eq. (60) can be separated to the angular and radial equations 
by letting 



* = Z[r)£[Q)e^ m ^ wt) , 
the separated angular equation can be found as 



(63) 



£" + cot 6£' 



X - 



(m + p cos i 
sin 2 # 



£ = 0, 



(64) 



where p = qQ and A is a separation constant. (From now on, a prime denotes the derivative 
with respect to its argument.) After setting the eigenvalue A = 1(1 + 1) — p 2 in Eq. (64), 
one can see that solutions to the angular part, £{6)e im ' f ', are the spin- weighted spheroidal 
harmonics p Yi m (9,ip) with spin- weight p 13]. 

On the other hand, before proceeding to the radial equation, one may rewrite the metric 
function f(r) in Eq. (28) as 



f(r) = -(r - r 2 )(r - n), 
r 



(65) 



where r 2 and r± denote the outer and inner horizons of the LDBHs, respectively. In the new 
form of the metric function Eq. (65), the physical parameters are 
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(66) 



r 2 = 7^ + Vc 2 - Aab^j , 
ri = ^ (c - Vc 2 - Aab^j , 



in which 



c = M = 4M and a = 



Ai 
A 2 A 2 ' 



(67) 



Since the algorithm in the calculations of the SCRSM cover only the outer region of the 
black hole (r > r 2 ), we must impose a condition in order to keep /(r) positive i.e. 6 > 0. 
Henceforth, one can derive the following radial equation as 



b(r - r 2 )(r - n)Z" + b(2r - r 2 - n)Z' + ( ^ 



r 2 a; 2 



(68) 



6(r — r 2 )(r — ri) A 2 

The above equation can be solved in terms of hypergeometric functions. Here, we give the 
final result as 



Z(r) = Ci(r - r 2 )^ r2 (r - n)-^ ri F 
C 2 {r - r 2 )- lQr2 {r - nY'^F 



a, 6; c; 



r 2 — r 
r 2 - ri 



+ 



a — c+1,6 — c+l;2 — c; 



r 2 — r 
r 2 - ri 



The parameters of the hypergeometric functions are 



(69) 



a = - + z(— + a), b = - + i(y - crj, and c = 1 + 2?wr 2 , 

where 



(70) 



Here, a is assumed to have real values. Furthermore, setting 



(71) 



x 



r-r 2 = exp(— ), 

W2 



(72) 
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one gets the behavior of the partial wave near the outer horizon (r — > r 2 ) as 

^ ~ C x e lw(x - () + C 2 e~ iui{x ~ t) . (73) 

One may infer the constants C\ and C 2 as being the amplitudes of the near-horizon outgoing 
and ingoing waves, respectively. 

In the literature, there exists a useful feature of the hypergeometric functions, which is 
a transformation of the hypergeometric functions of any argument (say z) to the hyperge- 
ometric functions of its inverse argument (l/z). The relevant transformation is given by 

Q 



F(a, b- c; z) = ^^-A(- Z )--- F (a, a + 1 - c; a + 1 - 6; l/z) (74) 
T{b)F(c — a) 

+ r^rr-'S ^)" 8 ^ & + 1 - S; & + 1 - a; l/z). 

The above transformation leads us to obtain the asymptotic behavior of the partial wave, 
easily. After applying the transformation to the general solution (69), we obtain the partial 
wave near-infinity as follows 



( r — r . \ —i&ri 

* ~ ^ — ; <{ B x exp % 



x 



- — (a + ujfjri) — tut 
7]r 2 



B 2 exp i 



x 



(— a + ufjri) — tut 

r]r 2 



(75) 



On the other hand, since we consider the case of r — > oo, the overall-factor term 



[r — r\ 

whence the partial wave (75) reduces to 



— ILUTl rsj 



exp i( 



xur\ . 
r 2 



(76) 



1 



yjr - r 2 



| Bi exp i 


X 


+ B 2 exp i 


X 


a — ut 


a — ut 








t]r 2 



(77) 



where Bi and B 2 correspond to the amplitudes of the asymptotic outgoing and ingoing 
waves, respectively. One can derive the relations between Bi,B 2 and C%, C 2 as follows 



£i _ Ci r ^) r (£-fe) | ^ T(2 -o)T(a-b) 



T(a)T(c-b) r(a-c+l)r(l-6) 



(78) 
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T(c)T(b-a) r(2-c)r(6 

+ U2" 



B — C 

1 r(b)r(c-a) "r(6-c + i)r(i-a) 

Hawking radiation can be considered as the inverse process of scattering by the black hole 
such that the outgoing mode at the spatial infinity should be absent [8j. Briefly B\ = 0, 
and it naturally yields the coefficient for reflection by the black hole as 



R 



\Ci\ 
IC,| 



r(c-6) |r(a)| 



r(i-fo) 



(79) 



\T(a-c+ 1) 



which is equivalent to 



R 

Thus the resulting radiation spectrum is 

N=(e% - l) 



COsh7T 




r 2 +ri \ 
r 2 —n ' 


cosh7r (a — 


COsh 7T 




fr 2 +n \ 
^r 2 —ri > 


coshvr(a+ |) 



(80) 



-l 



R 



Mi)' 



and finally, one can read the more precise value of the temperature as 



(81) 



T = w/ln 



COsh7T 


[ + w^ra+n 
o v r 2 — n 


cosh7r((T + f) 


COsh7T 


|"<7 -(^+n) 


cosh 7r (cr — ^) 



(82) 



This must be considered as the equilibrium temperature of the quantum field at the vacuum 
state valid for all frequencies. In the limit of ultrahigh frequencies (a ~ j), Eq. (82) reduces 
to 

T hiQh ~ limT ~ — r U r . -= ~ — , (83) 
' «>i ln[exp(^ xl 



which smears out the dependence and results in a pure thermal spectrum. One can 
immediately observe that Eq. (82) is independent from the horizons of the non-extreme 
LDBHs similar to the other 4-dimensional LDBH solutions possessing one horizon. 

But, contrary to the others js, 9], the resulting high frequency temperature Thigh Eq. (83) 
differs from the standard Hawking temperature Tjj [15J, which is computed as usual by 
dividing the surface gravity by 2n : 



T, 



H 



K 

2^ 



L 

Air 



Air 



(84) 
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Let us note that this same result for the Th can be obtained from Hawking's period argument 
of the Euclideanized line element. For this purpose we complexify time in (29) by t — > ir 
and rearrange the terms so that the line element reads in the form R 2 x S 2 , given by 

ds 2 ~ (£\ 2 + (pdrf + f 2 oh (d9 2 + sin 2 9dtf) . (85) 

Here r 2 h stands for the value of the radial coordinate on the outer horizon (when it exists) 
and the constant S Q reads 



4Ai 



\iA 2 



^M + M 2 -^^-2v) 



(86) 



which relates to the period of the angle r, upon the overall multiplication by S . Since T H 
is the inverse of the period we obtain 

1 



T, 



H 



2vr 



S 



(87) 



which is identical with (84), valid for double- horizon Hawking temperature. In order to find 
the vacuum 'in' and 'out' states for the scalar field we have to choose the metric such that 
the surface gravity and mass of the black hole both vanish. This can be done from (66), by 
choosing c 2 = 4ab first, to make an extremal LDBH (with zero temperature), and next, to 
let c — > to make the mass also zero. These conditions cast our LDBH metric into 



ds 2 = -brdt 2 + ^- + A 2 rdn 2 . 
or 

By simple arrangement this vacuum metric transforms into 

rfs 2 = p 2 (_ dT 2 + dx 2 + dQ 2j 

where 

r = e Px , t = ^-r, p = Aei x , (3 = AVb. 
b 

The massless Klein-Gordon equation V 2 $ = 0, with $ = ^ takes the form 



(88) 



(89) 



1 

? 



d TT -d xx + ^ + £(£ + l))V = 0. 



The vacuum 'in' and 'out' solutions for the scalar field are 

1 



-i(/3crx+uit) 

Vr- 
1 



(90) 



(91) 



out ~ —^e. 



i(f3ox—ujt) 
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where a has the meaning from (71). Once these states propagate from vacuum they turn 
into thermal states as described above. 

So the question arises here as in which case does the temperature Eq. (82) matches 
with the value of Tjj in Eq. (84)? The answer is absolutely related with the value of 
physical parameter, a . Let us assume that the value of the parameter a is so great that it 
predominates term ) (but ^( ^rj^q ) is still comparable with a) in the expression of 

the temperature (82). Unless this assumption is not violated, the corresponding limit of T 
will be Th- In summary, 



exp 7T 


<T-\ 


b » T2 


+■'■1 

— ri 


ex P 7r(a + f) 


exp 7i 


a — 


Ul /T-2 4 
b V T2- 


-n)" 


exp 7T (cr - |) 



Tff ~ lim T ~ lo/ In < ± ^ > (92) 

[*>f<i*£>]>i 

~ ^ ~ A ( i _ !1) 

In {exp 2vr [f (a±a)] exp 2tt (f ) } 4?r 

Another question may immediately come out: how does a maintain its predomination 
against ^( ^^ )? To clarify the question, one can check Eq. (71) in order to see that the 
predomination of a strictly depends on negative values of X. However, this is possible only 
with the case of p 2 = q 2 Q 2 > 1(1 + 1). Hence, a significant remark is revealed that obtaining 
Th of the non-extreme LDBHs from the SCRSM, the only possibility is to consider charged 
scalar waves instead of chargeless ones. 

Furthermore, we want to serve most intriguing figures about the spectrum temperature 
Eq. (82). To this end, first we plot T versus frequency u of non-extreme LDBHs with 
r i, r 2 7^ for low and high |p|-values, and display all graphs in Fig. 1. As it can be seen 
from Fig. 1 in the high frequencies the thermal behaviors of the LDBHs with different \p\- 
values exhibits similar behaviors in which their temperatures approach to Thigh while u — > oo. 
The plot with low |p|-value in Fig. 1 does not behave like the Hawking temperature. On 
the other hand, the other plot in Fig. 1, which has high |p|-value represents the Hawking 
temperature Tjj in the low frequencies (u > 0). Beside this, once the parameter a is lost its 
predomination against j (;rz^), the latter plot increases to reach the T high with increasing 
frequency as well. In the case of the non-extreme LDBHs with T\ = 0, there is no difference 
between Thigh and Th because of Eq. (84), and at the low |p|-values the temperature T 
exhibits similar behavior as in the case n , r 2 ^ 0, which is the well-known thermal character 

n 

in the EMD theory [8j. By the way, one should exclude u = during the plotting of the 
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temperature. Because it causes uncertainty for the temperature Eq. (82) and physically this 
case is not acceptable since we consider the propagation of scalar waves. Fig. 2 is about the 
graph of T versus frequency u of non-extreme LDBHs with r\ = in a high |p| -value. In 
this figure, it is illustrated that by increasing the frequency from + , the temperature first 
starts from a constant value, which is T# and then makes a peak (not much higher than 
Th), and then decreases back to Th while uj — > oo. Rousingly, one can observe that the 
behavior of the graph in Fig. 2 is very similar to the graph obtained from the well-known 
Planck radiation formula, see for instance 16] . Besides, both Fig. 1 and Fig. 2 show us 



that whenever high |p|-values are present, the frequency of the scalar wave needed to detect 
the temperature of the LDBHs as to be the Hawking temperature Th can either be very 
high (only for r\ — case, which is already known before |8j) or low. The latter information 
about the relationship between Th and low frequencies is completely new for us, and may 
play crucial role for the thermal detection of the LDBHs in the future. 

V. SINGULARITY ANALYSIS 



In section II, we present a solution in 4— dimensional static spherically symmetric EMD 
theory that incorporates two Liouville-type potential terms coupled with gravity together 
with magnetically charged dilatonic parameters. We have clarified that the solution possess 
a central singularity which is a characteristic feature for spherically symmetric systems. In 
the solutions that admit black holes this singularity is clothed by horizons. However, there 
are cases that this singularity is not hidden behind a horizons. In such cases the singularity 
is called a naked singularity. 

In classical general relativity, singularities are described as incomplete geodesies. This 
simply means that the evolution of timelike or null geodesies is not defined after a finite 
proper time. There is a general consensus that a removal of classical singularities is not 
important only for quantum gravity but also for other fundamental theories. In view of 
this consensus, we are aiming to analyze whether these classical naked singularities that 
occur in the general solution described in Eqs. (16)-(19) and in its linear dilaton limit 
given in Eqs. (27) and (28), turn out to be "strong" or "smoothed out" when probed with 



quantum test particles. Our analysis will be based on the pioneering work of Wald 12] 



which was developed by Horowitz and Marolf (HM)[10|]. HM, have proposed a criterion 
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to test the classical singularities with quantum test particles that obey the Klein-Gordon 
equation for static spacetimes having timelike singularities. The criterion of HM has been 



applied successfully for several spacetimes 17H19| within the context of quantum mechanical 



concepts. Among the others, HM have already analyzed the quantum singularity for the 
extreme case of the charged dilatonic black hole in the absence of Liouville-type potentials. 
They confirmed that for a specific interval of dilaton parameter, the singularity is quantum 
mechanically regular. The brief review of the criterion is as follows. 

A scalar quantum particle with mass m is described by the Klein-Gordon equation 
(V^V^ — m 2 ) ip — 0. This equation can be written by splitting the temporal and spatial por- 
tion as ~q£t = —Aif), such that the spatial operator A is defined by A = —\ffD l (a/7A) + 
fm 2 , where / = with ^ the timelike Killing field, while Di is the spatial covariant 

derivative defined on the static slice S. Then, the Klein-Gordon equation for a free rela- 
tivists particle satisfies = y/A~Eip, with the solution ip (t) = exp [it-sj Ae) if> (0) • If the 
extension of the operator A is not essentially self-adjoint, the future time evolution of the 
wave function is ambiguous. Then, HM criterion defines the spacetime quantum mechan- 
ically singular. However, if there is only one self-adjoint extension, the operator A is said 
to be essentially self-adjoint and the quantum evolution if) (t) is uniquely determined by the 
initial condition. According to the HM criterion, this spacetime is said to be quantum me- 
chanically regular. Consequently, a sufficient condition for the operator A to be essentially 



self-adjoint is to investigate the solutions satisfying the following equation ( see Ref. 20] for 
a detailed mathematical background), 

A$±i$ = 0. (93) 
This equation admits separable solution and hence the radial part becomes, 

«9 2 1 d(fR 2 ) d<j) + , m 2 4 

^ + 7i?^-^--^ -y 0± > = ' (94) 

in which I (I + 1) > is the eigenvalue of the Laplacian on the 2— sphere. The necessary 
condition for the operator A to be essentially self adjoint is that at least one of the solutions 
to this equation fails to be of finite norm when r — > 0. In summery, the self adjointness of the 
operator A, implies the well-posedness of the initial value problem. Therefore, the suitable 
norm \\(f)\\ for this case is the Sobolev norm which is used first time within this context by 
Ishibashi and Hosoya [20| defined by, 
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Uf = q j j ' R 2 r l \<P\ 2 dndr + \j ' R 2 f 

where q 2 is a positive constant and d\x is the volume element on the unit 2— sphere. The 
regularity of the central singularity at r = in quantum mechanical sense requires that the 
squared norm of the solutions of the Eq. (94) should be divergent for each 1(1 + 1) and each 
sign of imaginary term. The norm \\(j)\\ is divergent for I (/ + 1) > if it is for / = 0, so 
essential self-adjointness will be examined for I — (S — wave) case. This implies essential 
self adjointness for the operator A. Furthermore, we assume, a massless case (i.e. m — 0), 
and ignoring the term ±i-4 ( since it is negligible near the origin). 
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dr 



dfidr, 



(95) 



A. A more general case: 

The general solution for any value of a 2 which is related to the dilaton parameters 71 
and 72 is given in the Eq. (18). Since this solution is complicated enough for integrability, 
we consider the specific values of a 2 = 3 and V\ = 0. Hence, the general solution becomes; 

f(r) = (r - r 2 ) (r - n) , (96) 

where 



M ± VM 2 - A ai a 2 

ri < 2 = 2a~ 2 ' (97) 

Q 2 Ai 1 fQ 2 X 2 T/ 

° 1 = W a2= 12[^- V2 



The extreme case occurs when M 2 = Aa\a 2 . In this case there is one horizon only and it 

M 
2a 2 - 



is given by r h = If M > \/M 2 — Aaia 2 and M 2 > Aaia 2 , this particular case admits 



two horizons given by f i j2 . However, if M 2 — Aa\a 2 < 0, no black hole forms and hence, the 
singularity at r = becomes naked. 

As a requirement of the HM criterion, the singularity at r = must have a timelike 
character. This can be checked if one introduces tortoise coordinate defined by r* = J y 
and take its limit as r — > 0. We found that the limit is finite. Therefore, the singularity is 
timelike. The solution for Eq. (94) is 
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0(r) 



A- 



ln 



r — r 2 
r — f\ 



(98) 



16a 2 (^2 - ^ly 

The first and the second terms of the squared norm (95) is finite, when r — > 0. 

Consequently, the operator A is not essentially self-adjoint and therefore, the central 
singularity r = 0, remain quantum mechanically singular. 



B. The linear dilaton case: 



The metric function for the linear dilaton case can be written as, (from Eq. (65)) 



f(r) 



(r - r 2 ) (r - r x ) , 



where 



n,2 = 



M ± VM 2 - 4ab 



2b 



Ai 

A 2 A 2 ' 



The naked singularity occurs when M 2 — Aab < 0. The tortoise coordinate = f y is finite 
and indicating a timelike character at r = 0. The Penrose diagram of this particular case 
is shown in Fig. 3-a. The radial part of the separable Eq. (94) has solution for the linear 
dilaton case as, 



(f>{r) 



In 



r 2 



(99) 



bA 2 (r 2 — r\) r — T\ 
The first and the second term of the squared norm defined in Eq. (95) is finite. Therefore 
the spacetime is quantum mechanically singular. For the double- horizon case, M 2 — Aab > 0, 
which implies r± ^ r 2 ^ ^ ri, the timelike singularity at r = is not naked, and its Penrose 
diagram is depicted in Fig. 3-b. However, for a special case Ai = 0, the solution to Eq. (94) 
is 



em 



bM 



r -r h 



(100) 



in which = ^. The first term of the squared norm (95) is finite, whereas the second term 
behaves as, 
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~ (ln|r|)| r=0 ->oo. (101) 

Hence, under the condition Ai = 0, the central classical singularity becomes quantum me- 
chanically non-singular. When we have a single-horizon, with the choice Ai = 0, for example, 
the singularity r = 0, is shown in the Penrose diagram (Fig. 3-c). 



C. Near horizon behaviors 



In order to study the global behavior of our solution, at least for specific choices of 
parameters, and to be able to sketch the Penrose diagrams, we cast the metric into the form 
apt for near horizons. With the choice V\ — V<i — our metric function f(r) takes the form 

'"•7^(*-&+%ni)> <«*> 



in which 

, 1-a 2 



1 + a 2 ' 

and 



1 < 5 < 1 (103) 



A 2 = T^A 2 g 2 . (104) 
l — o 

Upon the choice of parameters involved, we can have, double, single or no-horizon cases. By 
a redefinition for time, our line element reads, in brief, 

ds 2 = A 2 ds 2 (105) 

in which 



d . 2 = _ {r r-Mr r +} ^ + _ _r_ _ _ ^ 2 + ^ 



(r — rJ) (r — r + ) j +2 | /•' 

(r — r_) (r — r + ) 

and 

r± = T+J [ J Si^ ' ' (107) 

We note that the global structure of ds 2 is same with ds 2 , and therefore we analyze ds 2 . The 
singularity structure of (106) can be seen from the Kretchmann scalar-K, which reads 

r -2(«5+3) 6 ^ ±1 

limfT-^ (108) 
r -° 1 r" 8 , 8 = +l 
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r -2(S+l) S ^ ±1 

lim K ~ ^ ^ (109) 

r ^°° I constant, 5 = — 1 

We concentrate ourselves now to the near horizon geometry by the following reparametriza- 

tion, with new coordinates (f , t) 

r_ = r , r + = r o + eb , r = r + er, t = (110) 

in which r and bo are constants and e — > 0, is a small parameter. We obtain, upon relabeling 
f = r and t — t 

d~s 2 = - r{r ~j> o) dt 2 + , r ° + \ dr 2 + rl +5 dQ 2 . (Ill) 
r i+6 r(r-b ) 

In Fig. 4 we plot the Penrose diagrams for the specific cases b = 0, and b > 0. The case 
b < doesn't differ from the case of b = 0, and as a matter of fact this particular case 
corresponds to the BR limit, which is known to correspond to the near extreme geometry of 
the RN black hole. The more standard BR is obtained from the present one by the inversion 
r — >■ -. 

r 



VI. CONCLUSION 



We have shown that dilaton field with Liouville's potential interpolates between RN 
black hole and non-black hole BR solution. The general solution for the metric function 
suggests that dilatonic presence induces significant changes in the solutions; for example, 
asymptotically flat black holes become non-asymptotically flat . It is shown, through radial 
linear perturbation, that dilaton can add instabilities to the otherwise stable RN black hole 
whereas BR remains stable. From the thermodynamic point of view also, by invoking specific 
heat the system can be tested against stability and phase transition. In the non-extreme 
LDBH case, which is a particular solution of our general solution the statistical and the 
standard Hawking temperatures are compared and plotted. It has been pointed out that 
with charged scalar waves and spin-weighted coupling the two results match for the case 
of double horizons. We recall that in the single horizon case, in spite of the existence of 
a linear dilaton such a discrepancy does not arise. It is remarkable that the spin-weighted 
spheroidal harmonics serve to convert the diverging temperature spectrum into a finite 
one. The presence of dilaton makes the spacetime highly singular at r = 0. Whether these 
singularities are quantum mechanically singular also, or not, we send a quantum test particle 
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and apply the criterion due to Horowitz and Marolf. We find that under certain choice of 
our parameters the naked singularities create an infinite repulsive quantum potential so that 
the particle feels a regular space time. 
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FIGURE CAPTIONS 

Figure 1: Temperature T as a function of oo for the non-extreme LDBHs in the case 
of r 1? r 2 7^ 0. The plots are governed by Eq. (82). Different line styles belong to different 
|p|-values: Dotted line corresponds to \p\ = 0.5 (as an example of low |p|-values) and solid 
line is for \p\ = 10 (as an example for high |p|-values). The physical parameters in Eq. (82) 
are chosen as follows: I — 1, b — 1, A — 1, r\ — 0.5 and r 2 = 1. 

Figure 2: Temperature T as a function of oo for the non-extreme LDBHs in the case of 
r 2 7^ and r\ = 0, and when p has a high value. The plot is governed by Eq. (82). The 
physical parameters in Eq. (82) are chosen as follows: \p\ = 10, / = 1, b = 1, A = 1 and 
r 2 = l. 

Figure 3-a: Penrose diagram for no- horizon case, M 2 — Aab < in which r = is a 
naked singularity. 

Figure 3-b: Penrose diagram of the LDBH with two distinct horizons r\ ^ r 2 , where 
r = is a timelike singularity. 

Figure 3-c: Penrose diagram of the LDBH with a single horizon at r = r h . singular 
nature of r = is not affected. 

Figure 4-a: Penrose diagram for the line element (111) with b = 0. There is no 
singularity and no horizon. By inverting coordinate r — > - we obtain the standard BR 
diagram. We note that the choice bo < is also similar to this case. 

Figure 4-b: For bo > there is a horizon and the Penrose diagram is as shown with 
singularities at the null boundaries. By inversion as in (3-a) we interchange r = (r — > oo) 
with r — > oo (r = 0). 
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